Chapter Eight 


Time Domain Controller Design 


8.1 Introduction 

In this chapter we study the problem of controller design such that the desired 
system specifications arc achieved. Controller design is performed in the time 
domain using the root locus technique. Controller design techniques in the 
frequency domain, based on Bode diagrams, will be presented in Chapter 9. 
In this book we emphasize controller design in the time domain for the following 
reasons: (a) with the help of MATLAB very accurate results can be obtained 
for both desired transient response parameters and steady state errors; (b) while 
designing a controller in time system stability will can be easily monitored since 
the root locus technique is producing information about the location of all of the 
system poles so that one is able to design control systems that have a specified 
relative degree (extent) of stability; (c) controller design using the root locus 
method is simpler than the corresponding one based on Bode diagrams', (d) root 
locus controller design techniques are equally applicable to both minimum and 
nonminimum phase systems, whereas the corresponding techniques based on 
Bode diagrams are very difficult to use (if applicable at all) for nonminimum 
phase systems (Kuo, 1995). 

Before we actually introduce the root locus techniques for dynamic controller 
design (Section 8.5), in Section 8.2 we consider a class of static controllers 
obtained through the pole placement technique based on full state feedback. The 
main puipose of this controller is to stabilize the closed-loop system, and it can 
sometimes be used to improve the system transient response. 
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In some cases it is possible to achieve the desired system performance by 
changing only the static gain K. In general, as K increases, the steady state 
errors decrease, but the maximum percent overshoot increases. However, very 
often a static controller is not sufficient and one is faced with the problem of 
designing dynamic controllers. 

In Section 8.3 we present common controllers used in linear system control 
design. Two main classes of these controllers are discussed: PI and phase-lag 
controllers that improve steady state errors, and PD and phase-lead controllers 
that improve the system transient response. Combinations of these controllers, 
which simultaneously improve both the system transient response and steady 
state errors, are also considered. 

A simple class of dynamic feedback controllers can be obtained by feeding 
back the derivative of the output variables. These controllers, known as the 
rate feedback controllers, are presented in Section 8.4. It is shown that a rate 
feedback controller increases the damping ratio of a second-order system while 
keeping the natural frequency unchanged so that both the response maximum 
percent overshoot and the settling time are reduced. 

Actual controller design in the time domain is studied in Section 8.5. Design 
algorithms (procedures) are outlined for several types of controllers introduced in 
Section 8.3. The impact of particular controllers on transient response parameters 
and steady state errors is examined. 

Several controller design case studies involving real physical systems are 
presented in Section 8.6. A the end of this chapter, in Section 8.7, a MATLAB 
laboratory experiment is formulated, in which students are exposed to the problem 
of controller design for real physical control systems. 

Chapter Objectives 

This chapter presents systematic procedures for time domain controller de¬ 
sign techniques based on the root locus method. Students will learn how to 
design different types of controllers such that the closed-loop control systems 
have the desired steady state errors and transient response parameters (maximum 
percent overshoot, settling time). In addition, the eigenvalue (pole) placement 
technique, which for controllable systems allows location of system eigenvalues 
in any desired position in the complex plane, is fully explained for the case of 
single-input single-output systems. 
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8.2 State Feedback and Pole Placement 


Consider a linear dynamic system in the state space form 

x = Ax + Bu 

r (8-D 

y = Cx 

In some cases one is able to achieve the goal (e.g. stabilizing the system or 
improving its transient response) by using the full state feedback, which represents 
a linear - combination of the state variables, that is 

u = -Fx (8.2) 


so that the closed-loop system, given by 

i = (A - BF)x 
y = Cx 


(8.3) 


has the desired specifications. 

The main role of state feedback control is to stabilize a given system so 
that all closed-loop eigenvalues are placed in the left half of the complex plane. 
The following theorem gives a condition under which is possible to place system 
poles in the desired locations. 

Theorem 8.1 Assuming that the pair (A,B) is controllable, there exists a 
feedback matrix F such that the closed-loop system eigenvalues can be placed 
in arbitrary locations. 

This important theorem will be proved (justified) for single-input single¬ 
output systems. For the general treatment of the pole placement problem for 
multi-input multi-output systems, which is much more complicated, the reader 
is referred to Chen (1984). 

If the pair (A, b) is controllable, the original system can be transformed into 
the phase variable canonical form, i.e. it exists a nonsingular - transformation 


such that 
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where cq’s arc coefficients of the characteristic polynomial of A, that is 

det(AI — A) — s n o, n _\S n ^ -|- d n — 2 ^ n T ■ • ■ T ctis T cto (8-6) 
For single-input single-output systems the state feedback is given by 


u( z) = J\ Z1 - f 2 Z 2 - - fnZ-n - -f c Z (8.7) 


After closing the feedback loop with u( z), as given by (8.7), we get from (8.5) 


0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

(«0 + fl) 

-(«! +fz) 

(«> + f:s) ... 

(^n—1 4" fn) - 


If the desired closed-loop eigenvalues are specified by A'/. A.j ..... A;;(, then the 
desired characteristic polynomial will be given by 

(A-*{)(*-AS)-(*-*) <8 . 9) 

= * n + ati^ 1 + <C 2** 2 + ■ • ■ + 4s + og 


Since the last row in (8.8) contains coefficients of the characteristic polynomial 
of the original system after the feedback is applied, it follows from (8.8) and 
(8.9) that the required feedback gains must satisfy 


a 0 + fl — a 0 fl — a 0 — Go 

«1 + h - <4 => h - <4 - «1 


1 fn — tt-a — 1 ^ fn — ^n—1 ^n — 1 


The pole placement procedure using the state feedback for a system which is 
already in phase variable canonical form is demonstrated in the next example. 

Example 8.1: Consider the following system given in phase variable canon¬ 
ical form 
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It is required to find coefficients /i, / 2 , h such that the closed-loop system has 
the eigenvalues located at Af 2 = — 1 ± jl, A3 = — 5. The desired characteristic 
polynomial is obtained from (8.9) as 

A d (A) = (A + 5)(A + 1 + jl)(A + 1 — jl) = A 3 + 7A 2 + 12A + 10 

so that from (8.10) we have 

/1 = «o - «o = 10 - 2 = 8 

h = 4 - «r = 12 - 5 = 7 

h = 4 ~ a-2 = 7 - 10 = -3 

o 

In general, in order to be able to apply this technique to all controllable 
single-input single-output systems we need to find a nonsingular transformation 
which transfers the original system into phase variable canonical form. This 
transformation can be obtained by using the linearly independent columns of the 
system controllability matrix 

C = b:Ab:A 2 b: - iA ?i 3 b 

It can be shown (Chen, 1984) that the required transformation is given by 

Q = [qi q 2 ••• q n ] (8.11) 

where 

q™ = b 

q n _i = Aq n + a n _ iq n = Ab + a n _ ib 

q^-2 — A.q^_^ T u ri _ 2 qri — A'b -)- Ab T u ri _ 2 b (8.12) 

qi = Aq 2 + aiq„ = A n_1 b + a n _\ A 7i_2 b H-f «ib 

where a -’s arc coefficients of the characteristic polynomial of matrix A. After 
the feedback gain has been found for phase variable canonical form, f c , in the 
original coordinates it is obtained as (similarity transformation) 


f = fcQ 


1 


(8.13) 
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Example 8.2: Consider the following linear system given by 
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The characteristic polynomial of this system is 

det(AJ - A) = s 3 + 1.3s 2 + 33s + 13 


Its phase variable canonical form can be obtained either from its transfer function 
(see Section 3.1.2) or by using the nonsingular (similarity) transformation (8.11). 
The system transfer function is given by 


c(sl — A) 1 b 


46s + 13 

s 3 + 13s 2 + 33s + 13 


Using results from Section 3.1.2 we are able to write new matrices in phase 
variable canonical form representation as 



0 

1 
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'O' 

A c = 
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-13 
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_ 1 _ 


c c = [13 46 0] 


The same matrices could have been obtained by using the similarity transfor¬ 
mation with 


A c = Q 1 AQ, b c = Q 3 b, c c = cQ 
with Q obtained from (8.11) and (8.12) as 


51 
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(8.14) 


Assume that we intend to find the feedback gain for the original system such 
that its closed-loop eigenvalues arc located at - 1 ,- 2 , —3, then 


A d (X) = (A + 1)(A + 2)(A + 3) = A 3 + 6A 2 + 11A + 6 
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From equation (8.10) we get expressions for the feedback gains for the system 
in phase variable canonical form as 

fi = 4 - a 0 = 6 - 13 = -7 

f 2 = a(-a i = ll- 33 = -22 

h■= 4 ~ a 2 = 6 - 13 = -7 

In the original coordinates the feedback gain is obtained from (8.13) 
f = f c Q 1 =[0.8149 -1.0540 5.7069] 


Using this gain in order to close the state feedback around the system we get 


x = 


-1.8149 

-0.6298 

-0.1851 


3.0540 

-0.8920 

-0.0540 


-5.7069 

-7.4139 

-3.2931 


x 


It is easy to check by MATLAB that the eigenvalues of this systems are located 
at -1,-2, —3. 


o 

Comment: Exactly the same procedure as the one given in this section can 
be used for placing the observer poles in the desired locations. The observers 
have been considered in Section 5.6. Choosing the observer gain K such that the 
closed-loop observer matrix A — KC has the desired poles corresponds to the 
problem of choosing the feedback gain F such that the closed-loop system matrix 
A 7 — F 7 B 7 has the same poles. Thus, for the observer pole placement problem, 
matrix A should be replaced by A T , B replaced by C 7 and F replaced by K T . 
In addition, it is known from Chapter 5 that the observability of the pair (A, C) 
is equal to the controllability of the pair (A 7 '. C 7 ), and hence the controllability 
condition stated in Theorem 8.1—the pair (A,B) is controllable, which for 
observer pole placement requires that the pair (A r , C T ) be controllable—is 
satisfied by assuming that the pair (A,C) is observable. 


8.3 Common Dynamic Controllers 

Several common dynamic controllers appeal - very often in practice. They are 
known as PD, PI, PID, phase-lag, phase-lead, and phase-lag-lead controllers. In 
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this section we introduce their structures and indicate their main properties. In 
the follow-up sections procedures for designing these controllers by using the 
root locus technique such that the given systems have the desired specifications 
are presented. In the most cases these controllers are placed in the forward path 
at the front of the plant (system) as presented in Figure 8.1. 



Figure 8.1: A common controller-plant configuration 


8.3.1 PD Controller 

PD stands for a proportional and derivative controller. The output signal of this 
controller is equal to the sum of two signals: the signal obtained by multiplying 
the input signal by a constant gain K p and the signal obtained by differentiating 
and multiplying the input signal by Kd, he. its transfer function is given by 

G c (s) = K p + Kds (8.15) 

This controller is used to improve the system transient response. 

8.3.2 PI Controller 

Similarly to the PD controller, the PI controller produces as its output a weighted 
sum of the input signal and its integral. Its transfer function is 

G c (s) = K p + Ki- = KpS + Kl (8.16) 

s s 
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In practical applications the PI controller zero is placed very close to its pole 
located at the origin so that the angular contribution of this “dipole” to the root 
locus is almost zero. A PI controller is used to improve the system response steady 
state errors since it increases the control system type by one (see Definition 6.1). 


8.3.3 PID Controller 


The PID controller is a combination of PD and PI controllers; hence its transfer 
function is given by 


1 K; 


G c {s) — P-p T P d& T Pi — 


Ii. pS 


K d s 2 


(8.17) 


s s 

The PID controller can be used to improve both the system transient response and 
steady state errors. This controller is very popular for industrial applications. 


8.3.4 Phase-Lag Controller 


The phase-lag controller belongs to the same class as the PI controller. The 
phase-lag controller can be regarded as a generalization of the PI controller. It 
introduces a negative phase into the feedback loop, which justifies its name. It 
has a zero and pole with the pole being closer to the imaginary axis, that is 


Gc(s) 


f Pl\ S + Zi 
\ Zi ) s + Pr ’ 


z\ > pi > 0 


(8.18) 


arg G c (s) = arg(* L n ) arg (s + pi) = 6 Zl - 0 Pl < 0 

where p\fz\ is known as the lag ratio. The corresponding angles 0 Zl and 0 Pl 
arc given in Figure 8.2a. The phase-lag controller is used to improve steady state 
errors. 


8.3.5 Phase-Lead Controller 

The phase-lead controller is designed such that its phase contribution to the 
feedback loop is positive. It is represented by 

n , \ s J r z 2 n 

G c (s) = — -, p 2 > Z 2 > 0 

S + P2 


G c (s ) = arg (s + z 2 ) - arg (s + p>) = 0 Z2 - 0 P2 > 0 


(8.19) 
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where 0 Z2 and 0 P2 arc given in Figure 8.2b. This controller introduces a positive 
phase shift in the loop (phase lead). It is used to improve the system response 
transient behavior. 



Figure 8.2: Poles and zeros of phase-lag (a) and phase-lead (b) controllers 


8.3.6 Phase-Lag-Lead Controller 

The phase-lag-lead controller is obtained as a combination of phase-lead and 
phase-lag controllers. Its transfer function is given by 

G c (s) = p-2 > -2 > -1 > Pi > 0, Z\ Z'2 = P\P2 (8.20) 

(s+pi)(s + p 2 ) 

It has features of both phase-lag and phase-lead controllers, i.e. it can be used to 
improve simultaneously both the system transient response and steady state errors. 
Flowever, it is harder to design phase-lag-lead controllers than either phase-lag 
or phase-lead controllers. 

Note that all controllers presented in this section can be realized by using 
active networks composed of operational amplifiers (see, for example, Dorf, 1992; 
Nise, 1992; Kuo, 1995). 


8.4 Rate Feedback Control 

The controllers considered in the previous section have simple forms and in 
most cases they are placed in the forward loop in the front of the system to 
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be controlled. Another simple controller that is always used in the feedback 
loop is known as the rate feedback controller. The rate feedback controller is 
obtained by feeding back the derivative of the output of a second-order system 
(or a system which can be approximated by a second-order system, i.e. a system 
with dominant complex conjugate poles) according to the block diagram given 
in Figure 8.3. 

The rate feedback control helps to increase the system damping. This follows 
from the fact that the closed-loop transfer function for this configuration is given 
by 

= _ 4 _ = 4 K > n , 82n 

U(s) S 2 + 2(C + \K t Un)u n s + t o 2 s 2 T 2 ( c u n s + m 2 ’ 4 

Compared with the closed-loop transfer function of the second-order system 
without control (6.4), we see that the damping factor is now increased to 

(c = C + (8.22) 

Since the natural frequency is unchanged, this controller decreases the response 
settling time (see (6.20)). The system response maximum percent overshoot is 
also decreased (see Problem 8.3). 



Figure 8.3: Block diagram for a rate feedback controller 


Example 8.3: Design a rate feedback controller such that the damping ratio 
of the second-order system considered in Example 6.1 is increased to Q = 0.75 
and determine the system response settling time and maximum percent overshoot. 
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From Example 6.1 we know 

u> n = 2 rad/s, £ = 0.5, t s ~ 3 s, MPOS = 16.3% 

The gain for the rate controller is obtained from (8.22) as 

fO = —2(Cc - C) = ^2(0.75 - 0.5) = 0.25 
oj n 2 

The new values for the response settling time and maximum percent overshoot 
arc given by 

4c ~ ——— = ---= 2 s, MPOS = e i()()(%) = 2.83 % 

Qc^n 0.15 X 2 

It can be seen that both the system response settling time and maximum percent 
overshoot are reduced. 

o 


8.5 Compensator Design by the Root Locus Method 

Sometimes one is able to improve control system specifications by changing 
the static gain K only. It can be observed that as K increases, the steady 
state errors decrease (assuming system’s asymptotic stability), but the maximum 
percent overshoot increases. However, using large values for K may damage 
system stability. Even more, in most cases the desired operating points for the 
system dominant poles, which satisfy the transient response requirements, do not 
lie on the original root locus. Thus, in order to solve the transient response and 
steady state errors improvement problem, one has to design dynamic controllers, 
considered in Section 8.3, and put them in series with the plant (system) to be 
controlled (see Figure 8.1). 

In the following we present dynamic controller design techniques in three 
categories: improvement of steady state errors (PI and phase-lag controllers), 
improvement of system transient response (PD and phase-lead controllers), and 
improvement of both steady state errors and transient response (PID and phase- 
lag-lead controllers). Note that transient response specifications are obtained 
under the assumption that a given system has a pair of dominant complex 
conjugate closed-loop poles; hence this assumption has to be checked after a 
controller is added to the system. This can be easily done using the root locus 
technique. 
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8.5.1 Improvement of Steady State Errors 

It has been seen in Chapter 6 that the steady state errors can be improved by 
increasing the type of feedback control system, in other words, by adding a pole 
at the origin to the open-loop system transfer function. The simplest way to 
achieve this goal is to add in series with the system a PI controller as defined 
in (8.16), i.e. to get 


G c .(s)G(a) = G {s ) 

Since this controller also introduces a zero at —Ki/K p , the zero should be placed 
as close as possible to the pole. In that case the pole at p = 0 and the zero at 
2 m p act as a dipole, and so their mutual contribution to the root locus is almost 
negligible. Since the root locus is practically unchanged, the system transient 
response remains the same and the effect due to the PI controller is to increase 
the type of the control system by one, which produces improved steady state 
errors. The effect of a dipole on the system response is studied in the next 
example. 

Example 8.4: Consider the open-loop transfer functions 


G i{s) 


(f + 2 ) 

(s + l)(s + 3) 


and 




(s + 2)(s + 5) 

(•s + l)( s T 3)(s T 5.1) 


Note that the second transfer function has a dipole with a stable pole at —5.1. 
The corresponding step responses arc given in Figure 8.4. It can be seen from 
this figure that the system with a stable dipole and the system without a stable 
dipole have almost identical responses. These responses have been obtained by 
the following sequence of MATLAB instructions. 


numl=[1 2 ] ; 

denl=[1 4 3]; 
num2=[1 7 10]; 

dl=[1 1]; 
d2=[l 3]; 
d3=[1 5.1]; 
dl2=conv(dl,d2); 
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den2=conv(dl2,d3); 

[cnuml,cdenl]=feedback(numl,deni,1,1,-1); 

[cnum2,cden2]=feedback(num2,den2, 1,1,-1) ; 

t=0:0.1:2; 

step(cnuml,cdenl,t) 

step(cnum2,cden2,t) 



Figure 8.4: Step responses of a system without 
a stable dipole (a) and with a stable dipole (b) 

It is important to point out that in the case of an unstable dipole the effect 
of a dipole is completely different. Consider, for example, the open-loop transfer 
function given by 


r tsi = v A ; 

31 J (s + l)(s + 3)(s — 5.1) 

Its step response is presented in Figure 8.5b and compared with the corresponding 
step response after a dipole is eliminated (Figure 8.5a). In fact, the system without 
a dipole is stable and the system with a dipole is unstable; hence their responses 
are drastically different. Thus, we can conclude that it is not correct to cancel an 
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unstable dipole since it has a big impact on the system response. 



Figure 8.5: Step responses of a system without an 
unstable dipole (a) and with an unstable dipole (b) 


o 

Both the PI and phase-lag controller use this “stable dipole effect”. They do 
not change the system transient response, but they do have an important impact 
on the steady state errors. 

PI Controller Design 

As we have already indicated, the PI controller represents a stable dipole with a 
pole located at the origin and a stable zero placed near the pole. Its impact on the 
transient response is negligible since it introduces neither significant phase shift 
nor gain change (see root locus rules 9 and 10 in Table 7.1). Thus, the transient 
response parameters with the PI controller are almost the same as those for the 
original system, but the steady state errors are drastically improved due to the 
fact that the feedback control system type is increased by one. 
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The PI controller is represented, in general, by 

s+ 

G c (s) = K P —IQ « K P (8.23) 

where K p represents its static gain and Kij K p is a stable zero near - the origin. 
Very often it is implemented as 

G c (s) = (8.24) 


This implementation is sufficient to justify its main purpose. The design algorithm 
for this controller is extremely simple. 

Design Algorithm 8.1: 

1. Set the PI controller’s pole at the origin and locate its zero arbitrarily close 
to the pole, say 2 C = 0.1 or z c = 0.01. 

2. If necessary, adjust for the static loop gain to compensate for the case when 
Kj, is different from one. Hint: Use K p = 1, and avoid gain adjustment 
problem. 

Comment: Note that while drawing the root locus of a system with a PI 
controller (compensator), the stable open-loop zero of the compensator will attract 
the compensator’s pole located at the origin as the static gain increases from 0 to 
+ oo so that there is no danger that the closed-loop system may become unstable 
due to addition of a PI compensator (controller). 

The following example demonstrates the use of a PI controller in order to 
reduce the steady state errors. 

Example 8.5: Consider the following open-loop transfer function 

G(s) = _ — ^ + _ 

K ’ (s + 10)(s 2 + 2s + 2) 

Let the choice of the static gain K = 10 produce a pair of dominant poles 
on the root locus, which guarantees the desired transient specifications. The 
corresponding position constant and the steady state unit step error are given by 



= 0.25 
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Using a PI controller in the form of (8.24) with the zero at —0.1 (z c — 0.1), we 
obtain the improved values as K p — oo and e ss — 0. The step responses of the 
original system and the compensated system, now given by 


G c (s)G{s) 


100 + 0 . 1 ) 0 + 6 | 

s(s + 10)(s 2 + 2s + 2) 


are presented in Figure 8 . 6 . 



Figure 8.6: Step responses of the original (a) 
and compensated (b) systems for Example 8.5 

The closed-loop poles of the original system arc given by 

Ai = -9.5216, A 2j3 = -1.2392 ± j'2.6204 

For the compensated system they arc 

A lc = -9.5265, A 2c , 3c = -1.1986 ± j'2.6109 

Having obtained the closed-loop system poles, it is easy to check that the 
dominant system poles arc preserved for the compensated system and that the 




348 


TIME DOMAIN CONTROLLER DESIGN 


damping ratio and natural frequency are only slightly changed. Using information 
about the dominant system poles and relationships obtained from Figure 6.2, we 
get 

(u; n = 1.2392, Lo 2 n = ( 1.2392) 2 + (2.6204) 2 => cu 2 = 2.9019, C = 0.4270 
and 

( c u> nc = 1.1986, cu 2 c = (1.1986) 2 + (2.6109) 2 
> u 2 nc = 2.8901, Cc = 0.4147 

In Figure 8.7 we draw the step response of the compensated system over a 
long period of time in order to show that the steady state error of this system is 
theoretically and practically equal to zero. 

Figures 8.6 and 8.7 are obtained by using the same MATLAB functions as 
those used in Example 8.4. 

The root loci of the original and compensated systems are presented in 
Figures 8.8 and 8.9. It can be seen from these figures that the root loci arc 
almost identical, with the exception of a tiny dipole branch near the origin. 



Figure 8.7: Step response of the compensated system for Example 8.5 
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Phase-Lag Controller Design 

The phase-lag controller, in the context of root locus design methodology, is also 
implemented as a dipole that has no significant influence on the root locus, and 
thus on the transient response, but increases the steady state constants and reduces 
the corresponding steady state errors. Since it is implemented as a dipole, its 
zero and pole have to be placed very close to each other. 

The lag controller’s impact on the steady state errors can be obtained from the 
expressions for the corresponding steady state constants. Namely, from (6.31), 
(6.33), and (6.35) we know that 

K p = lim (if(s)G(s)}, K v = lim {sH(s)G(s)}, K a = lim {s 2 i7(s)G(s)} 

5—^0 5—^0 5—^0 k } 

and from (6.30), (6.32), and (6.34) we have 

1 ^ 1 2 
e SSstep - 1 + J(p : e SSramp - —, e ss parabolio - jt- 

For control systems of type zero, one, and two, respectively, the constants 
K P ,K V , and K a are all given by the same expression, that is 

Ki = K ——-, l — p,v,a (8.25) 

PiP2 ■■■ 

Consider, first, a phase-lag compensator of the form 

G c (s) = —^ z c > p c > 0 (8.26) 

s + p c 

If we put this controller in series with the system, the corresponding steady state 
constants of the compensated system will be given by 

Z\ Zo ■ ■ ■ Z c Z c , , 

Ri c = R - = Ri—, Ic. I •• p. c. u (8.27) 

PlP2 ■■'Pc Pc 

In order to increase these constants and reduce the steady state errors, the ratio 
of z c jp c should be as large as possible. Since at the same time ,rv : must be close 
to p c (they form a dipole), a large value for the ratio $ C /Pc can be achieved if 
both of them are placed close to zero. For example, the choice of = 0.1 
and p c = 0.01 increases the constants Ki, l = p. v. a. ten times and reduces the 
corresponding steady state errors ten times. 
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Now consider a phase-lag controller defined by (8.18), that is 


G c (4 


(Pc\s+Z c 

\z c ) s + p c ' 


z c > p c > 0 


This controller will change the value of the static gain K by a factor of pjz c , 
which will produce a movement of the desired operating point along the root 
locus in the direction of smaller static gains. Thus, the plant static gain has 
to be adjusted to a higher value in order to preserve the same operating point. 
The consequence of using a phase-lag controller as defined in (8.18) is that the 
same (desired) operating point is obtained with higher static gain. We already 
know that by increasing the static gain, the steady state errors are reduced. In 
this case, the static gain adjustment has to be done by choosing a new static 
gain K = Kz c /p c . Note that the effects of both phase-lag controllers (8.18) and 
(8.26) arc exactly the same, since the gain adjustment in the case of controller 
(8.18) in fact cancels its lag ratio p c /z c . 

The following simple algorithm is used for phase-lag controller design. 

Design Algorithm 8.2: 


1. Choose a point that has the desired transient specifications on the root locus 
branch with dominant system poles. Read from the root locus the value 
for the static gain K at the chosen point, and determine the corresponding 
steady state errors. 

2. Set both the phase-lag controller’s pole and zero near - the origin with the 
ratio zjpc obtained from (8.27) such that the desired steady state error 
requirement is satisfied. 

3. In the case of controller (8.18), adjust for the static loop gain, i.e. take a 
new static gain as A = Kz c jp c . 


The next example demonstrates the controller design procedure with a phase- 
lag compensator according to the steps outlined in Design Algorithm 8.2. 

Example 8.6: The steady state errors of the system considered in Example 
8.5 can be improved by using a phase-lag controller of the form 

- 0.1 


G c {s) = 


0.01 


Since z c /p c = 10, the position constant is increased ten times, that is 


K* c = A' — = 3 x 10 = 30 

Pc 
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so that the steady state error due to a unit step input is reduced to 


e 


SS step 


1 

1 Kp C 


— = 0.03226 
31 


It can be easily checked that the transient response is almost unchanged; in 
fact, the dominant system poles with this phase-lag compensator are —1.2026 ± 
j'2.6119, which is very close to the dominant poles of the original system (see 
Example 8.5). 


o 


Example 8.7: Consider the following open-loop transfer function 


G(s)H(s) 


K(s + 15) 

s(s + 20)(s 2 + 4s + 8 ) 


Let the choice of the static gain K — 20 produce a pair of dominant poles on 
the root locus that guarantees the desired transient specifications. The system 
closed-loop poles for K = 20 arc given by 

Ai j2 = -0.5327 ± j'2.2024, A 3 = -2.9194, A 4 = -20.0153 


so that for this value of the static gain K the dominant poles exist, i.e. the 
absolute value of the real part of the dominant poles (0.5327) is about six times 
smaller than the absolute value of the real part of the next pole (2.9194), which 
is in practice sufficient to guarantee poles’ dominance. Since we have a type one 
feedback control system, the steady state error due to a unit step is zero. The 
velocity constant and the steady state unit ramp error are obtained as 

20 X 15 15 1 ro 

K v = -= — => e ss = —— = 0.53 

20 X 8 8 r '‘ mr K v 

Using the phase-lag controller with a zero at —0.1 (z c = 0.1) and a pole at 
— 0.01 (j> c = 0.01), we get 

_ 2 C 150 

K vc — K v — „ e ssc rarnp — 0.053 

Pc 8 

It can be easily shown by using MATLAB that the ramp responses of the original 
and the compensated systems arc very close to each other. The same holds for 
the root loci. Note that even smaller steady state errors can be obtained if we 
increase the ratio z c /p c , e.g. to z c jp c = 100 . 


o 
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8.5.2 Improvement of Transient Response 

The transient response can be improved by using either the PD or phase-lead 
controllers. In the following, we consider these two controllers independently. 
However, both of them have the common feature of introducing a positive phase 
shift, and both of them can be implemented in a similar manner. 

PD Controller Design 

The PD controller is represented by 

G c (s) = s + z c , z c > 0 (8.28) 

which indicates that the compensated system open-loop transfer function will 
have one additional zero. The effect of this zero is to introduce a positive phase 
shift. The phase shift and position of the compensator’s zero can be determined 
by using simple geometry. That is, for the chosen dominant complex conjugate 
poles that produce the desired transient response we apply the root locus angle 
rule given in formula (7.10) and presented in Table 7.1 as rule number 10. This 
rule basically says that for a chosen point, on the root locus the difference 
of the sum of the angles between the point s,{ and the open-loop zeros, and the 
sum of the angles between the point s,j and the open-loop poles must be 180°. 
Applying the root locus angle rule to the compensated system, we get 

m n 

G c { s d)G{ s d ) — ^{$d + %) + ^{ s d + z i) — 4-{ s d + Pi) — 180° (8.29) 

i =1 i =1 


which implies 

m n 

^{ s d + z c ) — 180° — ^{sd + Zi) + ^-{sd. + Pi) = o c (8.30) 

i =1 i =1 

From the obtained angle Z(s# + z c ) the location of the compensator’s zero is 
obtained by playing simple geometry as demonstrated in Figure 8.10. Using this 
figure it can be easily shown that the value of z c is given by 

A - (Oand c -i V'l 8 2 ) (8-31) 

tana c V / 


An algorithm for the PD controller design can be formulated as follows. 
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Figure 8:10 Determination of a PD controller’s zero location 


Design Algorithm 8.3: 

1. Choose a pair of complex conjugate dominant poles in the complex plane 
that produces the desired transient response (damping ratio and natural 
frequency). Figure 6.2 helps to accomplish this goal. 

2. Find the required phase contribution of a PD regulator by using formula 
(8.30). 

3. Find the absolute value of a PD controller’s zero by using formula (8.31); 
see also Figure 8.10. 

4. Check that the compensated system has a pair of dominant complex conjugate 
closed-loop poles. 

Example 8.8: Let the design specifications be set such that the desired 
maximum percent overshoot is less than 20% and the 5%-settling time is 1.5 s. 
Then, the formula for the maximum percent overshoot given by (6.16) implies 


/ = In {05} => c 


' In 2 {05} 

7T 2 + In 2 {05} 


0.456 


We take £ = 0.46 so that the expected maximum percent overshoot is less than 
20%. In order to have the 5%-settling time of 1.5 s, the natural frequency should 
satisfy 

•3 3 

is ~ U n « — = 4.348 rad/s 

(u n £L 
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The desired dominant poles are given by 

Sd = A d = ~(u n ± j^nV 1 - C 2 = -2.00 ± j3.86 
Consider now the open-loop control system 

G( s ) - _ — + 1Q ) _ 

U (s + l)(s + 2)(s + 12) 

The root locus of this system is represented in Figure 8.11a. 



-50 -45 -40 -35 -30 -25 -20 -15 -10 -5 0 

Real Axis 


Figure 8.11: Root loci of the original (a) and compensated (b) systems 

It is obvious from the above figure that the desired dominant poles do not belong 
to the original root locus since the breakaway point is almost in the middle 
of the open-loop poles located at —1 and —2. In order to move the original 
root locus to the left such that it passes through Sd, we design a PD controller 
by following Design Algorithm 8.3. Step 1 has been already completed in the 
previous paragraph. Since we have determined the desired operating point, Sd, we 
now use formula (8.30) to determine the phase contribution of a PD controller. By 
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MATLAB function angle (or just using a calculator), we can find the following 
angles 

£{s d + z\) — 0.4495 rad, + pi) = 1.8243 rad 

Z(sd + P 2 ) = 1.5708 rad, + p 3 ) = 0.3684 rad 

Note that MATLAB function angle produces results in radians. Using formula 
(8.30), we get 

l{s d + z c ) = ir- 0.4495 + 1.8243 + 1.5708 + 0.3684 
= 0.1723 rad = 9.8734° = a 0 

Having obtained the angle a c , the formula (8.31) produces the location of the 
controller’s zero, i.e. rn = 24.1815, so that the required PD controller is given by 

G c (s) = s + 24.1815 

The root locus of the compensated system is presented in Figures 8.11b and 
8.12b. It can be seen from Figure 8.12 that the point s d = — 2 ± j3.86 lies on 
the root locus of the compensated system. 



Figure 8.12: Enlarged portion of the root loci in the neighborhood of the 
desired operating point of the original (a) and compensated (b) systems 
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At the desired point, Sd, the static gain K, obtained by applying the root 
locus rule number 9 from Table 7.1, is given by K — 0.825. This value can be 
obtained either by using a calculator or the MATLAB function abs as follows: 

dl=abs(sd+pl); 
d2=abs(sd+p2); 
d3=abs(sd+p3); 
d4=abs(sd+z1) ; 
d5=abs(sd+zc) ; 

K=(dl*d2 *d3)/(d4*d5) 

For this value of the static gain K, the steady state errors for the original and 
compensated systems are given by e ss = 0.7442, e ssc = 0.1074. Note that in 
the case when z c > 1, this controller can also improve the steady state errors. 
In addition, since the controller’s zero will attract one of the system poles for 
large values of K, it is not advisable to choose small values for z c since it may 
damage the transient response dominance by the pair of complex conjugate poles 
closest to the imaginary axis. 

The closed-loop step response for this value of the static gain is presented 
in Figure 8.13. It can be observed that both the maximum percent overshoot and 
the settling time are within the specified limits. 



Figure 8.13: Step response of the compensated system for Example 8.8 
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The values for the overshoot, peak time, and settling time arc obtained by 
the following MATLAB routine: 

[yc,xc,t]=step(cnumc, cdenc) ; 

% t is a time vector of length i=73; 

% cnumc = closed-loop compensated numerator 
% cdenc = closed-loop compensated denominator 
plot(t,yc); 

[ymax, imax]=max(yc); 

% ymax is the function maximum; 

% imax = time index where maximum occurs; 

tp=t(imax) 

essc=0.1074; 

yss=l-essc; 

os=ymax-yss 

% procedure for finding the settling time; 
delt5=0.05*yss; 
i = 7 3; 

while abs ( (yc(i)-yss))<delt5; 

i=i-i; 

end; 
ts=t (i) 

Using this program, we have found that t s — 1.125 s and MPOS — 20.68%. 
Our stalling assumptions have been based on a model of the second-order system. 
Since the second-order systems are only approximations for higher-order systems 
that have dominant poles, the obtained results are satisfactory. 

Finally, we have to check that the system response is dominated by a 
pair of complex conjugate poles. Finding the closed-loop eigenvalues we get 
— —11.8251, A 2j 3 = —2.000 ± j3. 8600, which indicates that the presented 
controller design results are correct since the transient response is dominated by 
the eigenvalues A 2 , 3 . 


o 


Phase-Lead Controller Design 

The phase-lead controller works on the same principle as the PD controller. It 
uses the argument rule, formula (7.10), of the root locus method, which indicates 
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the phase shift that needs to be introduced by the phase-lead controller such that 
the desired dominant poles (having the specified transient response characteristics) 
belong to the root locus. 

The general form of this controller is given by (8.19), that is 

, v S + Z C „ 

G c [s) — : ■ Pc ^ ~c 0 

S + Pc 

By choosing a point Sd for a dominant pole that has the required transient response 
specifications, the design of a phase-lead controller can be done in similar fashion 
to that of a PD controller. First, find the angle contributed by a controller such 
that the point Sd belongs the root locus, which can be obtained from 

^Gc(sd) = 180° - lG(s d ) (8.32) 


that is 

m n 

0 C = £{s d + z c ) — ^-{sd + Pc) = 180° — ^-{sd + z{) + £{sd + Pi) (8.33) 

i =1 i =1 

Second, find locations of controller’s pole and zero. This can be done in many 
ways as demonstrated in Figure 8.14. 



Figure 8.14: Possible locations for poles and zeros of 
phase-lead controllers that have the same angular contribution 

All these controllers introduce the same phase shift and have the same impact on 
the transient response. However, the impact on the steady state errors is different 
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since it depends on the ratio of z c /p c . Since this ratio for a phase-lead controller 
is less than one, we conclude from formula (8.27) that the corresponding steady 
state constant is reduced and the steady state error is increased. 

Note that if the location of a phase-lead controller zero is chosen, then simple 
geometry, similar to that used to derive formula (8.31), can be used to find the 
location of the controller’s pole. For example, let —z c 3 be the required zero, then 
using Figure 8.14 the pole —p c 3 is obtained as 

Pc3 - (u n + u n y/l- C 2 tan (0 C - ip + 7t/2) (8.34) 

where ip — Z(.s„ + z c 3 ). Note that ip > 0 C . 

An algorithm for the phase-lead controller design can be formulated as 
follows. 

Design Algorithm 8.4: 

1. Choose a pair of complex conjugate poles in the complex plane that produces 
the desired transient response (damping ratio and natural frequency). Figure 
6.2 helps to accomplish this goal. 

2. Find the required phase contribution of a phase-lead controller by using 
formula (8.33). 

3. Choose values for the controller’s pole and zero by placing them arbitrarily 
such that the controller will not damage the response dominance of a pair of 
complex conjugate poles. Some authors (e.g. Van de Verte, 1994) suggest 
placing the controller zero at —(LJ n . 

4. Find the controller’s pole by using formula (8.34). 

5. Check that the compensated system has a pair of dominant complex conjugate 
closed-loop poles. 

Example 8.9: Consider the following control system represented by its 
open-loop transfer function 

GTA — _ ^ _ 

[ J (> + 10)(s 2 + 2s + 2) 

It is desired that the closed-loop system have a settling time of 1.5 s and a 
maximum percent overshoot of less than 20%. From Example 8.8 we know that 
the system operating point should be at *,/ = —2 ± j3.86. A controller’s phase 
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contribution, obtained from formula (8.33) is 

9 C = 7T - 0.7676 + 0.4495 + 1.9072 + 1.7737 
= 6.5044 rad = 0.2213 rad = 12.6769° 

Let us locate a zero at —15 (z c — 15), then by (8.34) the compensator’s pole 
is at —p c — —59.2025. The root loci of the original and compensated systems 
are given in Figure 8.15, and the corresponding step responses in Figure 8.16. 



Figure 8.15: Root loci for the original (a) and compensated (b) systems 

It can be seen that the root locus indeed passes through the point — 2 ± y'3.86. 
For this operating point the static gain is obtained as K — 101.56; hence the 
steady state constants of the original and compensated systems are given by 
K p — 30.468 and K pc — K p (z c /p c ) = 7.7196, and the steady state errors are 
e ss = 0.0317, e ssc = 0.1147. Figure 8.16 reveals that for the compensated system 
both the maximum percent overshoot and settling time are reduced. However, 
the steady state unit step error is increased, as previously noted analytically. 

Consider now another phase-lead compensator with a zero set at —9. From 
(8.34) we get p c — 15.291. The root locus of the compensated system with a 
new controller is given in Figure 8.17. 
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The static gain at the desired operating point — 2 ± j 3.86 is K — 41.587, 
and hence the steady state errors arc e ss — 0.0742, e ssc — 0.11986. The 
step responses of the original and compensated systems, for K — 41.587, are 
presented in Figure 8.18. 



Figure 8.18: Step responses of the original (a) and compensated 
(b) systems with the second controller for Example 8.9 

It can be seen that this controller also reduces both the overshoot and settling 
time, while the steady state error is slightly increased. 

We can conclude that both controllers produce similar transient characteristics 
and similar steady state errors, but the second one is preferred since the smaller 
value for the static gain of the compensated system has to be used. The 
eigenvalues of the closed-loop system for K — 41.587 arc given by 

A ic = -12.4165, A 2c = -10.8725, A 2Cj3c = -2.000 ± j3.8600 

which indicates that the response of this system is still dominated by a pair of 
complex conjugate poles. 


o 
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Remark: In some applications for a chosen desired point, sj, the required 
phase increase, 0 C , may be very high. In such cases one can use a multiple 
phase-lead controller having the form 

G?ead(s) = ( '“I"') , Pc>Z c > 0 (8.35) 

\S+PcJ 

so that each single phase-lead controller has to introduce a phase increase of 0 c /n. 


8.5.3 PID and Phase-Lag-Lead Controller Designs 

It can be observed from the previous design algorithms that implementation of a 
PI (phase-lag) controller does not interfere with implementation of a PD (phase- 
lead) controller. Since these two groups of controllers arc used for different 
puiposes—one to improve the transient response and the other to improve the 
steady state errors—implementing them jointly and independently will take care 
of both controller design requirements. 

Consider first a PID controller. It is represented as 

K s 2 4- K/ $ _i_ jC 

Gpid(s) — K p + KdS H-- = Kd -—-— 

s s (8.36) 

= Kd(s + z Cl ) ~ ^ — Gpd(s)Gpi(s ) 

which indicates that the transfer function of a PID controller is the product of 
transfer functions of PD and PI controllers. Since in Design Algorithms 8.1 and 

8.3 there are no conflicting steps, the design algorithm for a PID controller is 
obtained by combining the design algorithms for PD and PI controllers. 

Design Algorithm 8.5: PID Controller 

1. Check the transient response and steady state characteristics of the original 
system. 

2. Design a PD controller to meet the transient response requirements. 

3. Design a PI controller to satisfy the steady state error requirements. 

4. Check that the compensated system has the desired specifications. 
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Example 8.10: Consider the problem of designing a PID controller for the 
open-loop control system studied in Example 8.8, that is 

Q(s) - _ _ + 10 ^ _ 

(s + l)(s + 2)( a + 12) 

In fact, in that example, we have designed a PD controller of the form 

G pd (s) = s + 24.1815 


such that the transient response has the desired specifications. Now we add a 
PI controller in order to reduce the steady state error. The corresponding steady 
state error of the PD compensated system in Example 8.8 is e ssc = 0.1074. 
Since a PI controller is a dipole that has its pole at the origin, we propose the 
following PI controller 

s + 0.1 



In comparison to (8.36), we arc in fact using a PID controller with K,j = 1, 
z cl =24.1815, z C 2 =0.1. The corresponding root locus of this system compen¬ 
sated by a PID controller is represented in Figure 8.19. 



Figure 8.19: Root locus for the system from 
Example 8.8 compensated by the PID controller 
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It can be seen that the PI controller does not affect the root locus, and hence 
Figures 8.11b and 8.19 are almost identical except for a dipole branch. 

On the other hand, the step responses of the system compensated by the 
PD controller and by the PID controller (see Figures 8.13 and 8.20) differ in 
the steady state parts. In Figure 8.13 the steady state step response tends to 
y ss = 0.8926, and the response from Figure 8.20 tends to 1 since due to the 
presence of an open-loop pole at the origin, the steady state error is reduced 
to zero. Thus, we can conclude that the transient response is the same one as 
that obtained by the PD controller in Example 8.8, but the steady state error is 
improved due to the presence of the PI controller. 



Figure 8.20: Step response of the system from 
Example 8.8 compensated by the PID controller 


o 

Similarly to the PID controller, the design for the phase-lag-lead controller 
combines Design Algorithms 8.2 and 8.4. Looking at the expression for a phase- 
lag-lead controller given in formula (8.20), it is easy to conclude that 

Glag/leadi^s) — O7 l tad ('^ ) 


(8.37) 
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The phase-lag-lead controller design can be implemented by the following al¬ 
gorithm. 

Design Algorithm 8.6: Phase-Lag-Lead Controller 


1. Check the transient response and steady state characteristics of the original 
system. 

2. Design a phase-lead controller to meet the transient response requirements. 

3. Design a phase-lag controller to satisfy the steady state error requirements. 

4. Check that the compensated system has the desired specifications. 

Example 8.11: In this example we design a phase-lag-lead controller for a 
control system from Example 8.9, that is 

G(a) = 7 - A ^ + 6) -- 

v ; (s + 10)(s 2 + 2s + 2) 

such that both the system transient response and steady state errors arc improved. 
We have seen in Example 8.9 that a phase-lead controller of the form 

r ta — 5 + 9 

itad[ ! s + 15.291 

improves the transient response to the desired one. Now we add in series with 
the phase-lead controller another phase-lag controller, which is in fact a dipole 
near the origin. For this example we use the following phase-lag controller 

s + 0.1 
s + 0.01 



so that the compensated system becomes 


G(s) = G(s)G c (s) 


K(s + 6) (s + 9) (a+ 0.1) 

(a + 10)(s 2 + 2s + 2) (a + 15.291) (s + 0.01) 


The corresponding root locus of the compensated system and its closed-loop step 
response are represented in Figures 8.21 and 8.22. We can see that the addition of 
the phase-lag controller does not change the transient response, i.e. the root loci 
in Figures 8.17 and 8.21 are almost identical. However, the phase-lag controller 
reduces the steady state error from e ss j ea d = 0.11986 to c ss „ a(; /,.= 0.01344 
since the position constant is increased to 

0.1 _ 41.587 x 9 x 0.1 _ ^ Aon 
^pjag/lead ~ - 1Q x 2 X 15 X 0.01 ~ 
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so that 


1 


ss,lag/lead l _L f r 

-L i p la 


= 0.01344 


p, lag/lead 


O 


8.6 MATLAB Case Studies 

In this section we consider the compensator design for two real control systems: 
a PD controller designed to stabilize a ship, and a PID controller used to improve 
the transient response and steady state errors of a voltage regulator control system. 

8.6.1 Ship Stabilization by a PD Controller 

Consider a ship positioning control system defined in the state space form in 
Problem 7.5. The open-loop transfer function of this control system is 

= 0-8424 

s{s + 0.0546)(s + 1.55) 

The root locus of the original system is presented in Figure 8.23a. 



Figure 8.23: Root loci for a ship positioning control 
problem: (a) original system, (b) compensated system 
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It can be seen that this system is unstable even for very small values of the 
static gain. Thus, the system transient response blows up very quickly due to the 
system’s instability. Our goal is to design a PD controller in order to stabilize 
the system and improve its transient response. Let the desired operating point be 
located at sj — —0.2 ±j 0.3, which implies w n — 0.3606 rad/s and ( — 0.5547. 
We find from (8.30) that the required phase shift is a c — 72.0768°, and from 
(8.31) the location of the compensator zero is obtained at —0.297. Thus, the PD 
compensator sought is of the form 

G c (s) = s + 0.297 

It can be seen from Figure 8.23 that the root locus of the compensated system 
indeed passes through the point Sd — —0.2 ± j 0.3 and that the compensated 
system is stable for all values of the static gain. The static gain at the desired 
operating point is given by K Sd — 0.6258 and the corresponding closed-loop 
eigenvalues at this operating point are Ai c = —1.2046, A 2 C , 3 C = — 0.2± j'0.3. In 
Figure 8.24 the unit step response of the compensated system is presented. 



Figure 8.24: Step response of a ship positioning compensated control system 

It is found that y max — 1.2863, t p — 7.3043 s, and t s — 12.7826 s. From 
the same figure we observe that the steady state error for this system is zero. 
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which also follows from the fact that the system open-loop transfer function has 
one pole at the origin. 


8.6.2 PID Controller for a Voltage Regulator Control System 

The mathematical model of a voltage regulator control system (Kokotovic, 1972) 
is given in Section 6.7. The open-loop transfer function of this system is 


G( S ) 


154280 

(s + 0.2)(s + 0.5)(s + 10)(s + 14.28)0 + 25 ) 


The corresponding root locus is presented in Figure 8.25. Since one of the 
branches goes quite quickly into the instability region, our design goal is to move 
this branch to the left so that it passes through the operating point selected as 
= —1± jl. For this operating point, we have w n = \/2 rad/s and £ = 0.7071 
so that the expected maximum percent overshoot and the 5%-settling time of the 
compensated system are MPOS = 4.3214%, t s = 3 s. In addition, the design 
objective is to reduce the steady state error due to a unit step to zero. 



Figure 8.25: Root locus for a voltage regulator system 


We use a PID controller to solve the controller design problem defined above. 
The required phase improvement for the selected operating point is found by using 
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(8.30) as a c — 1.3658 rad = 78.2573°. From formula (8.31) the location of the 
compensator’s zero is obtained as — z c — —1.2079, so that the PD part of a PID 
compensator is 

G pd {s ) = s + 1.2079 

The branches of the root loci in the neighborhood of the desired operating point 
of the original and PD compensated systems arc presented in Figure 8.26. It 
can be seen that the compensated root locus indeed passes through the point 
Sd - -1 ± j 1. 



Figure 8:26: Root loci of the original (a) and PD (b) compensated systems 

The closed-loop unit step response of the system compensated by the PD 
controller is represented in Figure 8.27. Using the MATLAB programs given in 
Example 8.8, gives MPOS — 6.08%, t p — 2.1s, and t s — 3.5 s, which is quite 
satisfactory. However, the steady state unit step error is e ss p£> — 0.0808. Note 
that the static gain at the operating point, obtained by applying the root locus 
rule number 9 from Table 7.1, is K Sd — 4060.8. The closed-loop eigenvalues 
at the operating point are 

X lpD - -23.7027, A 2pd = -18.1675, A 3f3D 
A 4 i 5 pd = -0.997 ± jl.0011 


= -6.1105 




TIME DOMAIN CONTROLLER DESIGN 


373 


which indicates that the system has preserved a pair of dominant complex 
conjugate poles. 

In order to reduce this steady state error to zero we use a PI controller of 
the form 


s + 0.1 

s 

Since the compensated system open-loop transfer function now has a pole at the 
origin, we conclude that the steady state error is reduced to zero, which can also 
be observed from Figure 8.27. 




Figure 8.27: Step responses of PD and PID compensated systems 

The transient response specifications for the system compensated by the 
proposed PID controller are MPOS — 11.277%, t p — 2.1s, and t s — 3.1s. 
Thus, the proposed PI controller has slightly worsened the transient response 
characteristics. It is left to students, in the form of a MATLAB laboratory 
experiment, to check that the transient response specifications of the compensated 
system obtained by using PI controllers that have zeros located at —0.01 and 
— 0.001 are improved. 
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8.7 Comments on Discrete-Time Controller Design 

Similarly to the continuous-time controller design, the root locus method can 
be used for design of controllers (compensators) in the discrete-time domain. In 
Section 7.4 we have indicated that the discrete-time root locus method is identical 
to its continuous-time counterpart. Thus, the results presented in this chapter can 
be easily extended and used for the controller design of discrete-time systems. 
For more details, the reader is referred to Kuo (1992) for theoretical aspects, and 
to Shahian and Hassul (1993) for MATLAB discrete-time controller design. 


8.8 MATLAB Laboratory Experiment 

Part 1. Consider the control system given in Example 8.8, that is 

G( s ) - _ A (s + 10) _ 

1 J (s + l){s + 2)(s + 12) 

with the transient response requirements determined by a desired operating point 
located at sj = — 2 + j 3. 

(a) Design a PD controller such that the transient response requirements arc met. 

(b) Design a PID controller such that the steady state error is reduced to zero. 

(c) Design a phase-lead controller for the design problem defined in (a). 

(d) Add a phase-lag controller in series with the controller obtained in (c) in 
order to reduce the steady state error by 50%. 

Plot the root loci and unit step responses for (b) and (d) of both the original 
and compensated systems and compare the results obtained with PID and phase- 
lag-lead controllers. 

Part 2. Use MATLAB to design a phase-lag-lead controller for the control 
system of Example 8.9 such that t s < 2 s, MPOS < 10%, and e ssc < 0.01. 

Use the MATLAB program from Example 8.8 to find the actual response 
overshoot and settling time. Find the eigenvalues of the compensated system and 
check whether the system response is dominated by a pair of complex conjugate 
poles. 

Part 3. For the voltage regulator system considered in Section 8.6.2, use the 
same PD controller as in Section 8.6.2, but take for a PI controller the following 
forms: 
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(a) G PI {s ) = (s + 0.013/s. 

(b) G PI (s ) = (s + 0.001)/s. 

For both cases, find the unit step responses of the PID compensated systems 
and determine the transient response parameters. Compare the results obtained 
with the results from Section 8.6.2. 
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8.10 Problems 

8.1 Consider a single-input single-output system whose open-loop transfer func¬ 
tion is given by 


10 

s 4 + 3s 3 + 4s 2 + 1 


G{») = 
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Find the phase variable canonical form for this system and design the full 
state feedback static controller such that the system closed-loop poles arc 
located at Ai^ = -1 ± j'2, A3 = -3,A4 = -10. 

8.2 Find the feedback gain f such that the system 


X = 

'-1 

0 

2 

-2 

3 ' 

1 

x + 

"2' 

0 


1 

-1 

-5 


1 


with u — — fx has a pair of dominant complex conjugate closed-loop poles 
at — 1 ± j 2 and a closed-loop pole at —10. 

8.3 Plot the function 


and comment on the maximum percent overshoot dependence on the damp¬ 
ing ratio. 

8.4 Design a rate feedback controller such that for the second-order system 
represented by the closed-loop transfer function 

9 

s 2 + 3s + 9 

the maximum percent overshoot is 10%. Find the settling time of the 
compensated system. 

8.5 Consider the second-order system from Example 6.6. Its steady state error 
is e ss = 0.4. Design a PI controller to reduce the steady state error to zero. 
Use MATLAB to plot the root loci and unit step responses of the original 
and compensated systems. 

8.6 For the open-loop control system 


A (s + 

s(s + 20)(s 2 + 2s + 10) 

draw the root locus. Check that the static gain K — 3 produces a pair of 
complex conjugate dominant poles. Find the corresponding steady state 
errors and transient response parameters. Design a PI controller such 
that the steady state error is reduced to zero while the transient response 
characteristics are preserved. 
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8.7 Consider the controller design problem for the hydroturbine governors of a 
power system (Arnautovic and Skataric, 1991), represented by 



'-0.71 

0 

0 

0 

0 ' 


'0.71' 


0 

-2 

0 

0 

0 


2 

A = 

0.61 

1.28 

-1.46 

0.566 

0 

, B = 

0 


-0.18 

-0.37 

0.56 

-0.594 

-0.23 


0 


0 

0 

0 

314.16 

0 


0 


Assume that the output matrices are given by 

C = [0 0 1 0 1], D = 0 

Using MATLAB, perform the following: 

(a) Find its open-loop transfer function and the steady state error. 

(b) Suggest a phase-lag controller to reduce its steady state error ten times. 

(c) Locate the system operating point on the root locus approximately at 
—1,5 + jl and find the static gain at that point. 

8.8 Design a phase-lag controller for the system represented by 


w (s + l)(s + 5)(s + 10) 

which produces a steady state unit step error of less than 0.01. Take 
K = 200. 

8.9 Consider the synchronous machine from Section 7.5.2. Design a PI con¬ 
troller to reduce the steady state unit step error by 100%. Choose an op¬ 
erating point on the root locus and find the corresponding static gain. Use 
MATLAB to check that the closed-loop system is asymptotically stable. 

8.10 For the second-order system from Example 6.1, design a PD controller such 
that the compensated system has t s ~ 2 s and MPOS < 15%. 

8.11 Repeat Problem 8.10 with the following requirements t s ~ 2 s and 
MPOS < 10%. 

8.12 Solve Problem 8.11 using a phase-lead controller. 

8.13 Consider the problem of a PD controller design. Assume that a controller’s 

zero and desired operating point satisfy sj + z c ) > 90°. Derive a formula 

corresponding to (8.31). 
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8.14 Consider the phase-lead controller design for the system given in Example 
8.9. Use different values for the controller’s poles and zeros and examine 
their impact on the steady state errors. Suggest at least live different 
controllers. 

8.15 Derive formula (8.34) for finding the location of the pole of a phase-lead 
controller assuming that the location of its zero is chosen. 

8.16 Design a phase lead-controller such that 

G(s) = -- — ^ + 10) - 

w (s + 20)(s 2 + 2s + 2) 
has t s < 1.5 s and MPOS < 20%. 

8.17 Consider the voltage regulator control system from Section 8.6.2. Use a 
double PD controller of the form 

G pd (s) = (s + 1.2079) 2 
to compensate this system. 

(a) Draw the root locus of this system compensated with a double PD 
controller and compare it with the corresponding one from Section 
8.6.2 obtained using a single PD controller. 

(b) Choose the operating point for a pair of complex conjugate poles such 
that the damping ratio is the same as in Section 8.6.2, i.e. £ = 0.7071. 
Find the static gain at that point and the closed-loop eigenvalues. Does 
the compensated system preserve the response dominance of a pair of 
complex conjugate poles? 

8.18 Design a phase-lag-lead controller for the voltage regulator system from 
Section 8.6.2 such that MPOS < 5%, t s < 3 s, and e ss < 1%. 

8.19 Consider the ship position control system from Section 8.6.1. Design a 
phase-lead controller such that MPOS < 10% and t s < 10 s. 

8.20 Use MATLAB to find and plot the ramp responses of the original and 
compensated systems studied in Example 8.7. 

8.21 Consider the F-15 aircraft under supersonic flight conditions with its state 
space matrices given in Example 1.4. Note that this aircraft has one input 



TIME DOMAIN CONTROLLER DESIGN 


379 


and four outputs, where the outputs represent the state space variables. 
Thus, we arc able to get four transfer functions of the form 

Gi(s) — C,-(fiI - A) _1 B, i = 1,2,3,4 

(a) Using MATLAB, find all four transfer functions of the F-15 aircraft. 

(b) Plot the root loci for G,-(s), i = 1,2, 3,4. Comment on the stability 
properties of this aircraft with respect to the values of the static 
feedback gains Ki, i = 1,2, 3,4. 

(c) Find the closed-loop transfer functions with unit feedback and examine 
their stability. 

(d) Propose controllers which will stabilize all outputs of this aircraft. 

(e) For the proposed controllers that assure stabilization, find the closed- 
loop step responses and determine the response steady state errors 
and transient parameters. 

(f) If necessary, design dynamic controllers based on the root locus 
technique, as discussed in this chapter, such that steady state errors 
and transient response parameters arc improved. 

8.22 Repeat Problem 8.21 for the F-15 aircraft under subsonic flight conditions 
with the state space matrices given in Example 1.4. 1 


i 


Problems 8.21 and 8.22 can be assigned as either term papers or final projects. 



